The puzzle of 90° reorientation in the vortex lattice of borocarbide superconductors 
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We explain 90° reorientation in the vortex lattice of borocarbide superconductors on the basis of a 
phenomenological extension of the nonlocal London model that takes full account of the symmetry 
of the system. We propose microscopic mechanisms that could generate the correction terms and 
point out the important role of the superconducting gap anisotropy. 



Abrikosov vortices in type two superconductors repel 
each other and therefore tend to form two dimensional 
lattices when thermal fluctuations or disorder are not 
strong enough to destroy lateral correlations. In isotropic 
s-wave materials the lattices are triangular, however in 
anisotropic materials or for "unconventional" d-wave or 
p-wave pairing interactions less symmetric vortex lattices 
( VL) can_form as recent experiment on high T c cupratesu, 
SrRuOjn and borocarbides showed. The quality of sam- 
ples in the last kind of superconductors allows detailed 
reconstruction of the phase diagram by means of small 
angle neutron scattering, scanning tunnelling microscopy 
or Bitter decoration technique. For H\\c the presence of 
a whole series of structural transformations of VL was 
firmly established. At first, stable at high magnetic fields 
square lattice becomes rhombic, or " distasted triangu- 
lar", via a second order phase transitionou. Then, at 
lower fields, r-45° reorientation of VL relative to crystal 
axis occursQ'El. For H\\a a continuous lock- in phase tran- 
sition was predictedO. Above this transition apex angle of 
elementary rhombic cell of VL does not depend on mag- 
netic field, but below a critical field such a dependence 
appears. 

Theoretically the mixed state in nonmagnetic borocar- 
bide superconductors RN12B2C ', R = Y, Lu can be-,un- 
derstood in the frame of the extended London modelQ (in 
regions of the phase diagram close to H C 2 (T) -Line the ex- 
tended Ginzburg-Landau model can be useccm). So far 
this theory always provided qualitative and even quanti- 
tative description of phase transitions in VL and various 
other properties such as magnetization beha^ioiul, depen- 
dence of nonlocal properties on the disordert£l, etc. How- 
ever, recently another "reorientation" phase transition 
has been clearly observed in neutron scattering experi- 
ment on LuN 12B2C , which cannot be explained by the 
theory despite considerable efforts. When magnetic field 
of 0.3T was applied along the a axis of this tetragonal 
supecconductor sudden 90° reorientation of VL has been 
seertll. At this point a rhombic (nearly hexagonal, apex 
angle w 60°) lattice, oriented in such a way that the crys- 
tallographic axes are its symmetry axes, gets rotated by 
90°. Both initial and rotated lattices are found to coexist 
at the field range of about 0.1T wide around the transi- 
tion. Similar observations have been made in magnetic 
material ErNi2B 2 C. 

In this Letter we explain why the extended London 
model in its original form cannot generally explain even 



the existence of the 90° reorientation transition. The 
reason is that it possesses a "hidden" spurious fourfold 
symmetry preventing such a transition. Then we gener- 
alize the model to include the symmetry breaking effect 
and explain why the reorientation take place. Then we 
search for a microscopic origin of this effect. Using BCS 
type theory we find that anisotropy of the Fermi surface 
is ruled out due to smallness of its contribution. It is 
anisotropy of the pairing interaction that provides the 
required mechanism. We, therefore, suggest that there 
exist a correlation between the critical field of 90° reori- 
entation in VL and the value of the anisotropy of the 
gap. 

A convenieat starting point of any generalized "Lon- 
don" modelDlia is the linearized relation between the su- 
percurrent ji and the vector potential Aj : 



{A-K/c)jM = -KiM) A M)- 



(1) 



In the standard London limit the kernel Kij(q) is ap- 
proximated just by its q = limit, inverse mass matrix, 
while in the extended London model the quadratic terms 
of the expansion of the kernel near q = are also kepuj: 
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The significance of the quantity nij t ki is that it allows 
proper account of the symmetry of any crystal system 
while the first term does not guarantee this. At the 
same time it expresses nonlocal effects which are inher- 
ent to the electrodynamics of superconductors and below 
we call its component or their combination nonlocal pa- 
rameters. From its definition, ny^j = \ dq%q (ff ) 1 
is a tensor with respect to both the first and the second 
pairs of indices. However, the way riij^i transforms when 
the first and the second pairs of indices are interchanged 
is not obvious because the "origin" of these indices are 
quite different. The first pair (ij) comes, roughly speak- 
ing, from the variation of the free energy of the system " a 
superconductor in weakly inhomogeneous magnetic field" 
with respect to the vector potentia while the second pair 
(kl) comes from the expansion in vector q. Below we show 
that in general no symmetry n,j u = nu,ij is expected. 

The original derivation Eq. ||) from BCS theory in 
quasiclassical Eilenberger formulationQ produced a fully 
symmetric rank four tensor: Tij^jy ~ (viVjVkVi) with Vi 
being components of velocity of electrons at the Fermi 
surface. In this calculation independence of the gap 
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function on the orientation was assumed. Let us con- 
sider vortex lattice problem with this result. Specializing 
to tetragonal borocarbides, the number of independent 
component of tensor n ihkl is four: n aaaa< n aabb , n aacc 
and n cccc . In the case of external magnetic field oriented 
along a axis the free energy of VL, which is the relevant 
thermodynamic potential for a thin plate sample in per- 
pendicular external field, reads 

F=(B 2 /8ir)Y / ll + D(g x ,g y )r\ (3) 
D = \ 2 (m a gl + m c gl) + X 4 [n{m a gl + m c glf + dg 2 x gl] . 

Here B is magnetic induction and the summation runs 
over all vectors g of the reciprocal VL. The nonlocal 
parameters appearing in this equation have the form 
n = n aacc and d — n cccc m 2 + n aaaa vn} a 6n aacc m a Tn c . 
The free energy of Eq. (||) has been extensively studied 
numerically first minimizing it on the class of rhombic 
lattices with symmetry axes coinciding with the crystal- 
lographic axeaj and more recently by us for arbitrary 
lattices with one flux per unit cell. Despite the fact that 
great variety of vortex lattice transformation were iden- 
tified, no a 90° reorientation has been ever seen. The 
reason is quite simple: the considered free energy is ac- 
tually effectively fourfold symmetric. After rescaling the 
reciprocal lattice vectors 



9x — ► g x 



(4) 



the sum in Eq.(Q) becomes fourfold symmetric explic- 
itly. Based on this observation one concludes that ener- 
gies of the lattices participating in the 90° reorientation 
are equal exactly. Therefore no phase transition between 
them is possible in the framework of the extended Lon- 
don model of Eq.(|^) and further corrections are necessary 
to account for this transition. 

There might be a slight possibility that the observed 
90° reorientation presents the lock-in transition described 
in the beginning of this paper. For this to happen the 
rescaled square VL should looks almost hexagonal and, 
correspondingly, a particular value of masses asymmetry 
m a /m c = [cos (60°)/ cos (45 )] 2 = 1/2 is required. This 
is very different from the figures quoted in literatureO: 
m a /m c — 0.9/1.22 = 0.74. More importantly, according 
to this scenario one should see two degenerate lattices at 
small fields below the transition and only a single lattice 
at high fields above the transition which experimentally 
is clearly not the case. 

To explain 90° reorientation we proceed by correcting 
the model of Eq. (|3|). On general symmetry grounds 
for if 1 1 as one can expect more terms in the expression for 
D which describes vortex-vortex interactions. Given two 
fold symmetry of the present case we write down for D 
the expansion in Fourier series up to fourth harmonics, 
perform rescaling defined by Eq. (^) and obtain 

D eff = D (g) + D 4 (g) cos(4<^) + D 2 (g) cos(2<^), (5) 



where (f is the polar angle in the rescaled b— c plane. The 
quantity D from Eq. (|3j) produces only fourfold invariant 
terms: 



D (g) = X 2 g 2 + {n + d/8m a m c ) X 4 g 4 
Di{g) = - (d/8m a m c ) X 4 g 4 . 



(6) 
(7) 



The new term D 2 (g) expresses the effective fourfold sym- 
metry breaking. Experimentally, it should be small as 
indicated by recent success in qualitative understanding 
the angle dependence of magnetization of LuNi 2 B 2 Cu 
with field lying in the a — b plane on the basis of the 
theory without D 2 term. Accordingly, we can treat it 
perturbatively:F = + F (v " t) with 

F (0) - (B 2 /8tt) £ [1 + £> + £> 4 cos(Aip)}- 1 , (8) 



F {pert) = _ (B 2 /8n) 



D 2 COs(2lfi) 



[1 + D + D 4 cos(Aip)Y 



(9) 



where the summation is over g (see Eq. (Q)). The original 
degeneracy of two VL rotated by 90° with respect to each 
other is split now. To explain the 90° reorientation the 
sign of the perturbation should change at certain field 
B reo . Magnetic field influences the sum via constraint 
that area of the unit cell carries one fluxon. Roughly 
speaking D 2 (g) should change sign when g w yj B reo /§ . 
The simplest way to implement this idea is to write for 
D 2 (g) two lowest order terms in g : 



D 2 = W4g 4 



w 6 g 6 



(10) 



Quadratic term is not present since we have already 
rescaled it out in derivation of Eq. (H|). In principle 
the coefficient w§ can be derived from BCS similarly to 
n-ij^ki tensor within the framework of original extended 
London modelQ. Then it is proportional to the Fermi 
surface average of six components of Fermi velocity. To 
obtain IU4, however, the result riij ki ~ {viVjVkVi) of Ref. 
^ is not sufficient. Indeed, using general expression for 
nij,kl and repeating derivation of Eq. (^) from Eq. (0-|^) 
we see that 
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(11) 



In what follows we first demonstrate the presence of the 
first order phase transition in the model of Eq. ( |l0| ) and 
then provide a microscopical derivation of W4. 

The critical magnetic field of the 90° reorientation B reo 
depends only on the ratio r = —X 2 wq/w4 : . We deter- 
mined this dependence numerically using standard com- 
putational methods. At first, for a fixed B the equilib- 
rium form of VL unit cell was obtained by minimization 
of EaM. Then, the zero of the perturbation energy 
Eq. (g) was found. As usual] during the numerical cal- 
culations the cutoff factor exp (— £ 2 g 2 ) was introduced 
inside the above sums in order to properly account for 
the failure of the London approach in the vortex core. 
The calculated critical field is presented on Fig. 1 (we 
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used d = 0.05 and n = 0.015 typical for LuNi 2 B 2 C). We 
see that within the approximation of Eq. (|Io|) the 90° 
reorientation cannot happen at very low magnetic fields. 
For LuNi 2 B 2 C with A « 710 A, the field unit $ /(2ttA) 2 
is about 100 G. From the experimentally observed tran- 
sition field B reo = 2.95/cOaIil we estimate the relative 
strength of sixth and fourth order terms in D 2 (see Eq. 
©) as r = 0.036. 

To obtain wa we start by discussing a general pairing 
model which includes anisotropics in both the dispersion 
relation of electrons and the singlet pairing interaction 



H[ip] = hplie HV) - fi] ^ a + V, 



A 



(12) 



v = -- J i>t [i + s (iv)] i>l a ^- a [i + s (-»v)] i> a , 

where the summation over spin indices a =T,! is as- 
sumed. Here ip(x) an d ip^(x) are the electron destruction 
and creation operators, /i is chemical potential and A is a 
positive constant factorized from the pairing interaction 
for convenience. Dispersion relation e(k) and pairing in- 
teraction, of which <5(k) is a part, are usually defined in 
k-space. To treat magnetic field effects it is advanta- 
geous to define them in coordinate space. According to 
the rules of quantum mechanics, in the above functions 
of k we perform the replacements k — > — «V or k — > iV 
depending on whether derivatives act on tp or ip\ Then 
the standard minimal substitution — iV — ► II = — iV — A 
can be accomplished. This procedure, however, is not 
unique because the components of II do not commute 
with each other. Therefore, e(k) and S(k) are presented 
by their Taylor expansions and in those terms which con- 
tain mixed derivatives the symmetrization in IF is used. 

The kernel (q) from Eq. (jl]) is obtained by treating 
the effect of slowly varying magnetic field in terms of 
the linear response. The change in the Hamiltonian due 
to the presence of magnetic field Hi[ip,A] = H[ip,A] — 
H[ip, 0] is taken into account perturbatively. The result 
reads (see, for example, Ref. |l3|): 



d 2 H x 



dHi ggj 



, (13) 



where angular brackets denote the statistical average 
with unperturbed density operator exp (— H[ip, 0]/T) . 
Thus, we have to expand the functional Hi up to the 
terms quadratic in A. Because our aim is to calculate u>4 
we need only the coefficients of this expansion for A z and 
dAJdx. 

In its full generally the problem of Eq. (|lj) in mag- 
netic field is quite intractable and below we consider 
two particular cases which help us to estimate quan- 
titatively the magnitude of different contributions to 
W4: i) isotropic superconducting interaction, <5(k) = 0, 
and arbitrary dispersion relation e(k); ii) an example of 
weakly k dependent superconducting electronic interac- 
tion, <5(k) = Sgk 2 , and isotropic dispersion of the stan- 



dard form e(k) = k 2 /(2m). For simplicity in both cases 
a clean system was investigated. 
In the case i) we obtain 



Hi 



A z e >z -i(d x A z )^-(d 2 A z 



i>t [ A l e ,z 2 - iA z {d x A z ) e tZ 2 x 



-A z (d 2 A z 



£ ~2„2 



(d x A z ) 



2 £ Z 2 X 2 



Ipa, 



(14) 



where £ tZX means the second derivative of e(— iV) with 
respect to z and x components of the argument, and so 
on. The final results reads 



u>4 = 



4e 



2 dR 2 



(15) 



R = — cosh 
T 



2T 



E(k) = y/ (e(k) - /1) 2 + A 2 . 



At zero temperature R approaches zero exponentially 
and W4 vanishes. As temperature increases, W4 increases 
monotonically and reaches its maximal value at T = T c 
where it smoothly joins the corresponding component of 
q-dependent magnetic susceptibility tensor of the normal 
metal. For estimation we considered a simple dispersion 
relation e(k) = ^k 2 + jk z and assumed deviations from 
spherical Fermi surface to be small: a = am 2 fi -C 1. Ex- 
panding in a we obtain at T = T c that 
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2a%Jh^fi/2m. 



(16) 



where $0 = 2e/hc. This quantity is very small. In- 
deed, comparing it with the components of riij t kl pro- 
ducing contributions to Eq.(^|) we see that w^ s /n xxxx ~ 
a(A//i) 2 . Therefore in order to find an origin of 90° re- 
orientation one has to look elsewhere. 

The obvious possibility is to relax the assump- 
tion of the isotropic gap and turn to the case (ii). 
We calculated averages in Eq.(p"3|) using the 1/iV 
expansioro rather than the BCS approximation. The 
Hamiltonian Eq. (|l2|) becomes ip£ [e (— iV) — fj] i>% — 

[1 + 8 (iV)] tf*jt a [1 + 5 HV)] i> b a where N is 
number of (real or auxiliary) copies of the Fermi sur- 
face enumerated by a, b. The corresponding perturbation 
Hamiltonian found by the minimal substitution is 



A, 



1 

2m 



(d 2 z r a ) n] 



cc 



(17) 



Here the terms proportional to A 2 are omitted since they 
are local and cannot contribute to derivatives of K zz (q) 
with respect to q x required to obtain 104. For simpler sit- 
uations like the case (i) the leading order in 1/N expan- 
sion, with N set to 1, simply coincides with the BCS ap- 
proximation. The reason to resort to the 1/N expansion 
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is twofold. Firstly, the BCS expression for contains 
diagrams up to three loops (see Fig. 2c) which are very 
complicated. Secondly, unlike BCS, this nonperturbative 
scheme is systematically improvable. The last property is 
important when questions of principle are concerned. Af- 
ter observing that the order 1/N contributions, Fig. 2a, 
all vanish due to k ■<=>• — k asymmetry, we calculated the 
leading 1/N 2 contributions to the magnetic kernel, Fig. 
2b. At T = to leading order in Sq (further reducing 
number of integrals) the result reads 



5 8tt 
N 2 105 



(0 



4>o 



h 2 /j>/2m 



(18) 



where 5 = Sonifi is dimensionless gap anisotropy. There- 
fore in physical case of interest N — 1 we obtain 
w?4° p /fixxxx ~ S that is not necessary very small. This 
value is to be compared with w^ s originated from Fermi 
surface anisotropy which has huge suppression factor 
(/i/A) 2 . A noticeable angular dependence of the gap was 
indeed observed in the most recent Raxnan scattering ex- 
periments on Y and Lu borocarbideslij. 

To conclude, we found that the extended London 
model is incapable of explaining 90° reorientation in VL 
for H\\a because it produces an effective fourfold sym- 
metry of the free energy of VL. This symmetry becomes 
explicit after a rescaling transformation. We showed that 
in general case one should include into the extended Lon- 
don model correction terms for which Ji^jy ^ nu,ij (see 
Eq. (|J)). As a result, the true twofold symmetry of the 
system in magnetic field H\\a is restored and 90° reori- 
entation can be explained naturally. We demonstrated 
the two mechanisms that generate the correction terms: 
anisotropy of the Fermi surface and anisotropy of the 
superconducting gap, and showed that only the contri- 
bution of the latter one can lead to observable conse- 
quences. The investigation of vortex matter became re- 
cently a very sensitive tool to probe microscopic proper- 
ties of the superconductors. In this paper we employed 
it to infer qualitative and even quantitative information 
about pairing interaction by calculating nonlocal correc- 
tions to linear response. 

Note that inclusion of the correction terms will not 
change any conclusions of the extended London model for 
H\\c. On the other hand, nonzero "two fold symmetric" 
correction will lead jtp smearing, or even disappearance, 
of lock- in transition^ in VL for H\\a . Most probably 
it will be not possible to check this prediction in the 
same samples of LUN12B2C in which 90° reorientation 
was observed, because in this case the experimentally 
found opening angle of the unit cell of rhombic VLEJ 
indicates that the critical field of lock-in transition is far 
above H C 2- 
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FIG. 1. Critical field of 90° VL reorientation as a function 
of parameter r = —X 2 we/w4 (see Eq. (10)) 
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FIG. 2. a) Diagrams of the order 1/N, which vanish in 
our model, b) Diadrams of the order 1/N 2 contributing to 
W4. c) An example of complicated diagrams, which are of the 
order 1/V 3 and were neglected. 
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